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Abstract 

In this paper, we consider a multiple-input-multiple-output optical wireless communication (MIMO- 
OWC) system in the presence of log-normal fading. In this scenario, a general criterion for the design 
of full-diversity space code (FDSC) with the maximum likelihood (ML) detector is developed. This 
criterion reveals that in a high signal-to-noise ratio (SNR) regime, MIMO-OWC offers both large-scale 
diversity gain, governing the exponential decaying of the error curve, and small-scale diversity gain, 
producing traditional power-law decaying. Particularly for a two by two MIMO-OWC system with 
unipolar pulse amplitude modulation (PAM), a closed-form solution to the design problem of a linear 
FDSC optimizing both diversity gains is attained by taking advantage of the available properties on the 
successive terms of Farey sequences in number theory as well as by developing new properties on the 
disjoint intervals formed by the Farey sequence terms to attack the continuous and discrete variables 
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mixed max-min design problem. In fact, this specific design not only proves that a repetition code (RC) 
is the optimal linear FDSC optimizing both the diversity gains, but also uncovers a significant difference 
between MIMO radio frequency (RF) communications and MIMO-OWC that space dimension alone 
is sufficient for a full large-scale diversity achievement. Computer simulations demonstrate that FDSC 
substantially outperforms uncoded spatial multiplexing with the same total optical power and spectral 
efficiency, and the latter provides only the small-scale diversity gain. 

Index Terms 

Multiple-input-multiple-output (MIMO), optical wireless communications (OWC), log-normal fad¬ 
ing channels, linear space code, full diversity, repetition coding and maximum likelihood detector. 


1. Introduction 

In the past decade, the demand for capacity in cellular and wireless local area networks has 
grown in an explosive manner. This demand has triggered off an enormous expansion in radio 
frequency (RF) wireless communications. As an adjunct or alternative to RF communication, 
optical wireless communications (OWC), due to its potential for bandwidth-hungry applications, 
has become a very important area of research 0 -p^. The importance of OWC lies in the 
advantages of low cost, high security, freedom from spectral licensing issues etc. Furthermore, 
OWC links of practical interest involve satellites, deep-space probes, ground stations, unmanned 
aerial vehicles, high altitude platforms, aircraft, and other nomadic communication partners. 
Moreover, all these links can be used in both military and civilian contexts, or both indoor and 
outdoor scenarios in demand of high data rate. Therefore, OWC is considered to be the next 
frontier for net-centric connectivity for bandwidth, spectrum and security issues. 

However, some challenges remain, especially in the mobile or atmospheric environments. 
For high data rate OWC systems over mobile or atmospheric channels, robustness is a key 
consideration. In mobile links, there will be inevitable impairments such as terminal-sway, aerosol 


scattering and non-zero pointing errors [ 131-[ 161, etc. In addition, for atmospheric environments, 
atmospheric effects, such as rain, snow, fog and temperature variation, will affect the link 
performance. Therefore, in the design of OWC links, we need to consider these impairments- 
induced fading 0- This fading of the received intensity signal can be described by the log¬ 


normal (LN) statistical model [18|-|23|, which is considered in this paper. To combat fading, 
multi-input-multi-output (MIMO) OWC (MIMO-OWC) systems provide diverse replicas of 
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transmitted symbols to the receiver by using multiple receiver apertures with sufficient separation 
between each so that the fading for each receiver is independent of others. Such diversity can 
also be achieved by introducing a design for the transmitted symbols distributed over transmitting 
apertures (space) and (or) symbol periods (time). Full diversity is achieved when the decaying 
speed of the error curve for the coded MIMO-OWC system is maximized. 

Unfortunately, unlike MIMO techniques for radio frequency (MIMO-RF) communications 
with Rayleigh fading, there are two significant challenges in MIMO-OWC communications. 
The first is that there does not exist any available mathematical tool that could be directly 
applied to the analysis of the average pair-wise error probability (PEP) when LN is involved. 
In this scenario, it is indeed a challenge to extract a dominant term of the average PEP. Eet 
alone say how to achieve a full diversity gain. Here, it should be mentioned that there are 
really mathematical formulae in literature for numerically and accurately computing the integral 


involving EN [18|-[22|. However, it can not be used for the theoretic analysis on diversity. The 
second is a nonnegative constraint on the design of transmission for MIMO-OWC, which is a 
major difference between MIMO RE communications and MIMO-OWC. It is because of this 
constraint that the currently available well-developed MIMO techniques for RE communications 
can not be directly utilized for MIMO-OWC. Despite the fact that the nonnegative constraint can 
be satisfied by properly adding some direct-current components (DC) into transmitter designs so 


that the existing advanced MIMO techniques for RE communications such as orthogonal 
space-time block code (OSTBC) [|25|-p^ could be used in MIMO-OWC, the power loss arising 


from DC incurs the fact that these modified OSTBCs [34|-[36| in an EN fading optical channel 


have worse error performance than the RC [211, [37|-[39|. 

All the aforementioned factors greatly motivate us to develop a general criterion on the 
design of full-diversity transmission for MIMO-OWC. As an initial exploration, we consider 
the utilization of a space dimension alone, and intend to uncover some unique characteristics of 
MIMO-OWC. With this goal in mind, our main tasks in this paper are as follows. 

1) To establish a general criterion for the design of full-diversity space code (FDSC). To 
this end, our main idea here is that by fragmenting the integral region of the average 
PEP involving EN into two sub-domains adaptively with SNR, the dominant term will be 
extracted. With this, we will give a necessary and sufficient condition for a space code to 
assure full diversity. 
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2) To attain an optimal analytical solution to a specific two by two linear FDSC design 
problem. To do that, we formulate and simplify our optimization problem with continuous 
and discrete mixed design variables by taking advantage of available properties as well as 
by developing new properties on Farey sequences in number theory for our purpose. In 
fact, we will prove that the RC is the optimal linear FDSC for this specific scenario in the 
sense of the design criterion proposed in this paper. 

II. Channel Model And Space Code 

In this section, we first briefly review the channel model which is considered in this paper. 

Then, we propose the space coding structure and formulate the design problems to be solved. 


A. Channel Model with Space Code 


Let us consider an M x N MIMO-OWC system having M receiver apertures and N 
transmitter apertures transmitting the symbol vector s, {se},£ = 1, ■■■ , L, which are 

randomly, independently and equally likely, selected from a given constellation. To facilitate 
the transmission of these L symbols through the N transmitters in the one time slot (channel 
use), each symbol S£ is mapped by a space encoder to an x 1 space code vector (s^) and 
then summed together, resulting in an x 1 space codeword given by x = J2e=i where 

the n-th element of x represents the coded symbol to be transmitted from the n-th transmitter 
aperture. These coded symbols are then transmitted to the receivers through flat-fading path 
coefficients, which form the elements of the M x N channel matrix H. The received space-only 
symbol, denoted by the M x 1 vector y, can be written as 


y 



Hx -f n, 


( 1 ) 


where Pop is the average optical power of x and, the entries of channel matrix H are independent 
and LN distributed, i.e., hij = where Zij ~ Af [pij, f = 1, ■ ■ • , M, j = 1, ■ ■ ■ , N. 
The probability density function (PDF) of hij is 


fn (hij) = 


V^^TT hij (J ij 


exp 


(\]lhij 


(2) 


The PDF Of H is /h (H) = n£i 117=1 fn {K). 
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The signalling scheme of s is unipolar pulse amplitude modulation (PAM) to meet the 
unipolarity requirement of intensity modulation (IM), i.e., s G As an example, the 

constellation of unipolar 2^-ary PAM is B 2 P = {0,1, ■ ■ ■ , 2^ — 1}, where p is a positive 
integer. Then, the equivalent constellation of s is 5 = {s : Sj G B 2 P.,i = 1, ■■■ , T}, i.e., 
S = B^,. 

Furthermore, for noise vector n, the two primary sources at the receiver front end are due to 
noise from the receiver electronics and shot noise from the received DC photocurrent induced by 


background radiation [40|, |411. Although the signal intensity also results in shot noise, which is 
signal-dependent, the shot noise from hight-intensity background radiation dominates. Therefore, 
by the central limit theorem, this high-intensity shot noise for the lightwave-based OWC is 


closely approximated as additive, signal-independent, white, Gaussian noise (AWGN) [411 with 
zero mean and variance cr^. 

By rewriting the channel matrix as a vector and aligning the code-channel product to form a 
new channel vector, we have Hx = (Im 0 x^) vec (H), where ® denotes the Kronecker product 
operation and vec (H) = [hu, ..., / iiat , ..., Hmi, ■ • •, hMN]^- For discussion convenience, we 
call Im C)x^ a codeword matrix. Then, the correlation matrix of the corresponding error coding 


matrix between Im ® and 1m 0 x^ is given by 

(Im ® x"^ - 0 x^)^ ( 1 m ® x'^ - Im 0 x"^) = Ijw ® (ee^) 


(3) 


where e is the “distance” between distinct codewords x and x. All these non-zero e form an 
error set, denoted by S. 


B. Problem Statement 

To formally state our problem, we make the following assumptions throughout this paper. 


l^i=l 


= p, 


op‘ 


1) Power constraint. The average optical power is constrained, i.e., E 

Although limits are placed on both the average and peak optical power transmitted, in the 
case of most practical modulated optical sources, it is the average optical power constraint 
that dominates [ [421 . 

, since the noise variance per 


2) SNR definition. The optical SNR is defined by pop = 




dimension is assumed to be cr^. Thus, in expressions on error performance involved in the 
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squared Euclidean distance, the term p, in fact, is equal to 



(4) 


with optical power being normalized by Unless 


. Unless Stated otherwise, p is referred to as 


the squared optical SNR thereafter. 

3) Channel state information. Channel state information (CSI) at the receiver (CSIR) is known 
and CSI at the transmitter (CSIT) is unavailable. 

Under the above assumptions, our primary task in this paper is to establish a general criterion 
on the design of FDSC and solve the following problem. 

Problem 1: Design the space encoder F(-) subject to the total optical power such that 1) 
Vs G iS, F (s) meets the unipolarity requirement of IM; 2) Full diversity is enabled for the MF 
receiver. ■ 

Naturally, two questions come up immediately: 1) What is full diversity referred to as for 
MIMO-OWC ? 2) What is the design criterion ? So far, both questions remain open and thus, 
motive us to analyse the error performance in the ensuing section. 

III. Error Performance Analysis And Design Criterion for Space Code 

In this section, our purpose is to analyze the error performance of the MF detector and 
develop a design criterion for FDSC. To make our presentation more clear as well as our main 
idea more easily understandable, we begin with the symbol error probability (SEP) analysis on 
single-input-single-output OWC (SISO-OWC). 


A. SEP of SISO-OWC with PAM 



. In this specific 


case, it is known [431 that the average SEP for 2^-ary PAM is given by 



(5) 


where p denotes the squared optical SNR defined by Q, ci = 2 — 2^“^ and C 2 = 2 '^v-i • the 
sake of simplicity and without loss of generality, we assume p = 1 which is the case of on-off 
keying (OOK). 
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Lemma 1: For OOK, (p) is bounded by 

1 


p 2 \ 2 

In p —In — 

4ct2 


^ op 


^ < -Pe (p) 


=2 ^2 


, In P-In ^ 

< - (Inp) 


, (lip) 

+ O I e 


( 6 ) 


The proof of Lemma is postponed into Appendix 

In the proof of Lemma our selection of r is reasonable in finding the dominant term 
adaptively with SNR. Here, it should mentioned that the term “adaptive” is referred to the fact 
that the fragmentation of the integral region is done by a function of SNR. We also noticed that 
the splitting method is similar to those used in [|44|. However, the significant difference 
between our proposed technique and those used in [fT^, [|44| is that our proposed integral-splitting 


method is to extract the dominant behavior by selecting r adaptively with SNR. The SEP for 
SISO-OWC has been given and the adaptive fragmentation with SNR can help attack the integral 
involved in LN. In light of this, we extend this technique to the case of MIMO-OWC. 


B. PEP of MIMO-OWC 

This subsection aims at deriving the PEP of MIMO-OWC by means of the above-mentioned 
adaptive fragmentation with SNR and then, establishing a general design criterion for the space 
coded MIMO-OWC system. 


Given a channel realization H G 


s,MxN 


and a transmitted signal vector s, the probability of 


transmitting s and deciding in favor of s with the ME receiver is given by [451 

' d fe 


P(s ^ s|H) = Q 


where (e) = X]i=i (hfe) with hj = [ha, ■■■ , hiNf , i = 1, 
over H yields 


F is —)■ s') = 


J F(s^s|H)/h(H) dU. 


(7) 

M. Averaging 0 

( 8 ) 


Eor presentation clarity, we first give the definition of the dominant term. Suppose that there are 
two terms Ti{p) and T 2 {p), each of which goes to zero when p tends to infinity. We say that 
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Ti(p) is a dominant term in the sum of Ti{p) + T 2 (p) if limp_).oo T 2 (p)/Ti(p) = 0. 
the dominant term of we make an assumption for time being. Later on, we will 
this condition is actually necessary and sufficient for a full diversity achievement. 

Assumption 1: Any e G is unipolar without zero entry. 

We are now in a position to formally state the first main result in this paper. 
Theorem 1: Under Assumption P (s —>■ s) is asymptotically bounded by 


To extract 
prove that 


^ n \-MN Z.i=iZ-j=i 

-V- 

Pl(s-s>s) 


< P is —s') 


<a(lnp)-“^^e 


E m v-viv 

i=l z-i=l 




8^ 

'>‘3 


+0 e 


~ go .2 


(9) 


Pu{s->-s) 


where U = Ef=i 



( 471 )"^" exp (-^) 


,fc=i 



and 


Cu = 


(NPi) 


MN 


o TT^r T~[^ ! ' 

211i=ilL=i \T 


. Tl.(NP%Tl 
exp-In ' 


8 


M 


The detailed proof of Theorem is provided in Appendix 

With all the aforementioned preparations, we are able to give the general design criterion for 
FDSC of MIMO-OWC in the following subsection. 


C. General Design Criterion for FDSC 


The discussions in Subsection III-B tells us that Pu (s —s) is the dominant term of the upper- 


bound of P (s —s) in With this, we will provide a guideline on the space code design in 
this subsection. To define the performance parameters to be optimized, we rewrite Pu (s —s) 
as follows. 


Pu (s -)■ s) = CuGc (e) 


P 


In p 


if In 
4 


ivp2pn 

M 


-flnPaCe) 
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X. n \-MN I ^ 1 
X (In p) exp ( — — In 


P 


InV 


where V, (e) = and (e) = exp E*=i (In 


( 10 ) 

Afp2 n \ 


M 


Here, the following three factors dictate the minimization of Pu (s —)■ s): 


1) Large-scale diversity gain. The exponent f2 with respect to In governs the behavior 

of Pu ^ s). For this reason, is named as the large-scale diversity gain. The full 

large-scale diversity achievement is equivalent to the event that all the MN terms in 

^ ~ 'n!j=i offered by the N x M MIMO-OWC are fully utilized. On the other 

hand, O can be considered to be the reciprocal of the channel equivalent variance. In 

fact, if the channel is independently and identically distributed with variance being aj^, 
2 

then ^ In this case, the full large-scale diversity gain can be considered to be the 

reduction of the channel variance by MN, which is the maximum reduction amount. In the 
sense of this point, our definition of O being the large-scale diversity gain is parallel to the 
diversity order for STBC of MIMO-RF. Thus, when we design space code, full large-scale 
diversity must be assured in the first place to maximize the exponential decaying speed 
of the error curve. 


2) Small-scale diversity gain. Vg (e) = k: 


I Af 
^1 = 1 ^ 


is called small-scale diversity gain. 


which affects the power-law decaying in terms of miner’s (e) should be maximized 
to optimize the error performance of the worst error event. Since the small-scale diversity 
gain will affect the average PEP via the power-law decaying speed of the error curve, the 
small-scale diversity gain of the space code is what to be optimized in the second place 
to maximize the power-law decaying speed of the error curve. 

3) Coding gain. Gc (e) is defined as coding gain. On condition that both diversity gains 
are maximized, if there still exists freedom for further optimization of the coding gain, 
maxee£- Gc (e) should be minimized as the last step for the systematical design of space 
code. 

In what follows, we will give a sufficient and necessary condition on a full large-scale diversity 
achievement. From ( [T0| ), we know that Assumption is sufficient for FDSC. Here, from the error 
detection perspective, we prove that Assumptionis also necessary. Let us consider the following 
two possibilities. 

1) Bipolar error vector. Without loss of generality, assume there exists an error vector e 
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containing bipolar entries, i.e., e* < 0 and Cj > 0,z 7^ j,i,j G {1, ■ ■ ■ , N}. Then, there 
exists a positive vector ho = [ei, • • • , e*, ••• , cnY' such that (e^ho)^ = 0, where 
Ci = III addition, for any .^ > 0, = .^ho is positive and (e^h^)^ = 

(e^ho.^)^ = 0 = 0. In other words, there exist two distinct signal vectors x and x 

satisfying (x —x)^ho = 0. That is to say, this signal design is not able to provide the 
unique identification of the transmitted signals in the noise-free case, and, then the reliable 
detection of the signal will not be guaranteed, even in a sufficiently high SNR. Therefore, 
the maximal decaying speed of the corresponding error curve can not be achieved. As a 
matter of fact, in this case, full large-scale diversity gain can not be attained. 

2) Error vector with zero entries. If n entries of e are zero, with index being N—n+1, ■ ■ ■ , N 
without loss of generality, then, ■ Therefore, there 

must be Mn entries of the channel matrix H that makes no contribution to PER In other 
words, all the degree of freedoms offered by the N x M MIMO-OWC are not utilized 
and, as a consequence, full large-scale diversity can not be achieved. 

To sum up, if either e is bipolar or has zero-valued entries, the corresponding large-scale 
diversity gain will be less than MN. Hence, Assumption is sufficient and necessary for FDSC, 
which is summarized as the following theorem: 

Theorem 2: A space code enables full large-scale diversity if and only if Ve G e is unipolar 
without zero-valued entries or equivalently, Ve G S, ee^ is positive. ■ 

It is known that for MIMO-RF, the maximal decaying is achieved if and only if the difference 
matrix of any two distinct codeword matrix is full-rank [24|. However, from Theorem]^ we can 
see that a space-only transmission can assure the maximal decaying speed of the corresponding 
PEP, say, full large-scale diversity gain. The main reason for this significant difference is that the 
MIMO-OWC channels are nonnegative and the coding matrix is not necessary to be full-rank, 
which is verified by Theorem. 

Now, the questions raised at the end of Section |I^ is answered. With these results, we can 
proceed to design FDSC systematically in the following section. 


IV. Optimal Design of Specific Finear FDSC 

In this section, we will exemplify our established criterion in ( [T^ by designing a specific 
linear FDSC for 2 x 2 MIMO-OWC with unipolar pulse amplitude modulation (PAM). For this 
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particular design, a closed-form solution to the design problem of space codes optimizing both 
diversity gains will be obtained by taking advantage of the available properties on the successive 
terms of Farey sequences in number theory as well as by developing new properties on the 
disjoint intervals formed by Farey sequence terms. 


A. Design Problem Formulation 


Consider a 2 x 2 MIMO-OWC system with F (s) = Fs, where F = 


/ll /l2 
/ 2 I /22 


and 


T 

ee = 


eiCa 


. By Theorem 


6162 


ee^ should be positive to maximize the large-scale 


diversity gain. On the other hand, from me structure of ee^ and ( fTO] ), the small-scale diversity 
gain is Vg (e) = |eie 2 | under the assumption that CSIT is unknown. Therefore, to optimize the 
worst case over £, FDSC design is formulated as follows: 


max mm 6162 

/ll,/l2,/21,/22 e 


s.t. 


[61,62]^ G £Jij > 0,i, j e {1,2}, 


( 11 ) 


6162 > 0, /ll -f /12 -f /2I + /22 = 1- 

Our task is to analytically solve ( [TT] ). To do that, we first simplify ( [TT] ) by finding all the 
possible minimum terms. 


B. Equivalent Simplification of Design Problem 

For 2^-PAM, all the possible non-zero values of 6162 are 

6162 = (m/ll ± ^/i 2 ) (m /21 ± nf 22 ) fO,m,nE B 2 P. 


( 12 ) 


1) Preliminary simplification: After observations over ( [T^ , we have the following facts. 

1) Vm 7 ^ 0 ,m,n G B2P, it holds that 

(m/ll + nf 12) (m/21 + n/22) > /ii/2i. 

2) Vn 7^ 0 , m, 77 , G B2P, it is true that 

(m/ll + nfi2) (m/21 + nf22) > /i2/22- 
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3) V/c 7 ^ 0, 7 ^ 0, k,m,n G B 2 P, we have 

k (m/ll - n/ 12 ) (m /21 - n/ 22 ) 


> 1 . 


(m/ll - nfi2) (m/21 - n/22) 

So, all the possible minimum of 6162 in ca are /11/21, /12/22 and (m/n - n/12) (m/21 - ri/22), 
where — is irredueible, i.e., m ± n. These terms are denoted by 

-^10 = /12/22 X , Foi = /12/22, 


-^mn = / 12/22 ^m^ - 


m 


. /21 
/22 


n 


After putting aside the eommon term, / 12 / 22 , we can see that Fmn is the piecewise linear function 
of ^ and respectively. So, ( [TT] ) can be solved by fragmenting interval [0, 00 ) into disjoint 
subintervals. This fragmentation can be done by the breakpoints where F^n = 0. To characterize 
this sequence, there exists an elegant mathematical tool in number theory presented below. 

2) Farey sequences: First, we observe some specific examples of the breakpoint sequences. 
For OOK, the breakpoints j, cxo. For 4-PAM, they are Y,|,|,|,x,|,f,f,C) 0 . For 8-PAM, we 
have the breakpoint sequence with the former part being 

0111121231432534561 
1 ’ 7 ’ 6 ’ 5 ’ 4 ’ 7 ’ 3 ’ 5 ’ 7 ’ 2 ’ 7 ’ 5 ’ 3 ’ 7 ’ 4 ’ 5 ’ 6 ’ 7’1 

and the remaining being 

76547357275374567 
6 ’5’4’3’5’2’3’4’T’3’2’1’2’T’1’l’T’°° 

Through these special examples, we find that the series of breakpoints before 1/1 (such as the 

sequence in ( |13a| )) is the one which is called the Farey sequence [^. The Farey sequence for 

any positive integer k is the set of irreducible rational numbers | with 0 < a < b < k arranged 

in an increasing order. The series of breakpoints after f (such as the sequence in ( |13b[ )) is the 

reciprocal version of the Farey sequence. Thus, our focus is on the sequence before j. 


(13a) 


(13b) 


The Farey sequence has many interesting properties [46|, some of which closely relevant to 
our problem are given as follows. 

Proposition 1: If —, — and — are three successive terms of 5^, A; > 3 and tiL -c ^ 
then. 


1 ) min 2 — m 2 ni = 1 and mi + m 2 > k + 1. 


2 ) 


ni+n2 

mi+m2 


e ^ ) and ^ = ^^4^. 

V mi ' m3 ' ^ -\-m.o 


m 2 


mi+m3 * 
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However, having only Proposition is not enough to solve our design problem in ( [TT] ). We 
need to develop the other important new properties of Farey sequenees, whieh will be utilized 
in the FDSC design in the ensuing subseetions. 

Property 1: Given k > 3, assume — E 'Sk and If tix. and 

^ are successive, then, ^ and < ^. ■ 

m 2 mi+m3 — m 2 mo+m2 — mi 

The proof of Property is postponed into Appendix 
Property 2: Assume ^ e ^k-,k > 3 and ^ Then, 

^ m\ ’ m2 mi m2 


1 ) — < ^ 

mi mi+m2 


^ holds. 

m 2 


m2n2* 


> F 

mini m2n2* 


2) <^’' 

~ ^ ~ ml+mV Pmini “ Tm2n2- 

The proof of Property is given in Appendix 
Using Properties and we attain the following property. 

Property 3: If ^ and ^ are successive in and ^ then, F, 


F 


m2n2 


mi m 2 

are the two worst cases. 


mi ’ m 2 


and 


The proof of Property is provided in Appendix 


C. Techniques to Solve The Max-min Problem 

Thanks to Farey sequences, ( [TT] ) is transformed into a piecewise max-min problem with two 
objective functions. By solving this kind of problem, another main result of this paper can be 
formally presented as the following theorem. 

Theorem 3: The solution to ([TT]) is determined by 


FDSC — 


2 + 2P+1 



or 


FDSC 


2 + 2P+1 


2P 1 

2P 1 


(14a) 


(14b) 


The proof of Theorem [^ is postponed into Appendix [^ 
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Theorem uneovers the faet that the optimal linear spaee eoded symbols are aetually unipolar 
22p-ary PAM symbols, since B 22 P = {si + 2^S2 : si, S 2 ^ B 2 P }. Therefore, in fact, we have proved 


that RC |211 is optimal in the sense of the criterion established in this paper. 


V. Computer Simulations 

In this section, we carry out computer simulations to verify our theoretical results. We first 
examine the performance bounds given in Lemma 1 and Theorem 1. As shown by Figs, [^and 
it can be seen that in the high SNR regimes, the proposed upper-bounds have almost the 
same negative slope as those of simulated results. In other words, the proposed upper-bounds 
have captured the behavior of the error curve with respect to the decaying speeds. As shown by 
Figs. 0 and 1^ when SNR is sufficiently high, our proposed bounds have a horizontal shift to the 
right compared with the simulated results. Therefore, the tightness of the proposed asymptotical 
bounds is mainly dependent on the precise estimate of some constant independent of SNR. 



Fig. 1 

Error performance of SISO-OWC 


In the following, we simulate to verify our newly developed criterion in ( [T^ . In light of 
our work being initiative, the only space-only transmission scheme available in the literature 
is spatial multiplexing. Accordingly, we compare the performance of spatial multiplexing and 


FDSC specifically designed for 2 x 2 MIMO-OWC in Section IV In addition, we suppose that 
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1 bit pcu, Repetition code with M=1, N=2 



10 “ 


10 


15 


20 


25 


Squared Opticai SNR/dB 


Fig. 2 


Error performance of MIMO-OWC 


= 1,2 are independently and identically distributed and let an = = 0-21 = (^22 = cr. 

These schemes are as follows: 


1) FDSC. The optical power is normalized in such a way that Yl^ij=ifij = 2 yields 
^ Zli j=i fij^j = 1- From ( [T4l ), the coding matrix is | 



2) Spatial Multiplexing. We fix the modulation formats to be OOK and vary So the rate 
is 2 bits per channel use (pcu). The transmitted symbols si,S 2 are chosen from {0,1} 
equally likely. The average optical power is E [si + S 2 ] = 1- 
We can see that both schemes have the same spectrum efficiency, i.e., 2 bits pcu and the same 
optical power. Through numerical results, we have following observations. 

1) Substantial enhancement from FDSC is achieved, as shown in Fig. [s} For = 0.01, the 
improvement is almost 16 dB at the target bit error rate (BER) of 10“^. For cr^ = 0.5, 
the improvement is almost 6 dB at the target BER of 10“^. Note that the small-scale gain 
also governs the negative slope of error curve. The decaying speed of the error curve of 
FDSC is exponential in terms of In and that of spatial multiplexing is power-law with 
respect to p. The reason for this difference is that FDSCs are full-diversity guaranteed by 
the positive constraints in ( [TT] ), whereas spatial multiplexing does not satisfy the positive 
requirement in Theorem 
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2bit5 pcu, 2x 2 MIMO-OWC 



Fig. 3 

BER COMPARISONS OF FDSC AND SPATIAL MULTIPLEXING. 



Squared Optical SNR/dB 


Fig. 4 

BER COMPARISONS OF SPATIAL MULTIPLEXING. 


2) Spatial multiplexing presents only small-seale diversity gain illustrated in Fig.|^ By varying 
the varianee of H, we find that in the high SNR regimes, the error eurve deeays as p~^ 
as long as the SNR is high enough. From cr^ = 0.001 to cr^ = 0.1, the error eurve has 
a horizonal shift, whieh is the typieal style of MIMO RF p^. The reason is given as 
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I el 6162 \ 

follows. The equivalent spaee eoding matrix is e = , 61,62 G {0,±1} 

y 6162 el j 

with 6 ^ + e| 7 ^ 0. It should be noted that there exists two typieal error events: 6162 = —1 
and 6162 = 0. From the necessity proof of Theorem]^ for 6162 = —1, the attained large- 
scale diversity gain is zero, and at the same time, if 6162 = 0 with ef -f- 7 ^ 0 , then, 

the attained large-scale diversity gain is only two for 2 x 2 MIMO-OWC. Therefore, the 
overall large-scale diversity gain of spatial multiplexing is zero with small-scale diversity 
gain being attained. 

3) From Fig. and Q we notice that when a increases, the error performances of repetition 
codes will worsen. For this phenomenon, the reason is that for repetition codes, the 
large-scale diversity is given by = Ylij=i which dominates the exponential 
decaying speed. However, for spatial multiplexing, increasing cx^ will improves the error 
performance by providing a horizontal shift to the left. It is known that since an exact 
error probability formula for OWC over log-normal fading channels is indeed hard to be 
obtained particularly for MIMO-OWC, it is very challenging to theoretically prove that 
the relationship between the error performance of spatial multiplexing and Intuitively 
speaking, the relationship between the error performance of spatial multiplexing and cx^ 
follows the radio frequency MIMO. That is to say, is the average electrical power of 
In h. That’s why increasing will improve the error performance of spatial multiplexing, 
as shown by Figs. and 

We can see that the performance gain of MIMO-OWC if any will become larger against 
increasing SNR in a high enough regime. This implies that a slight improvement in the 
coding structure will result in a significant enhancement of error performance in the high SNR 
regimes instead of only an horizonal shift to the left. Unique characteristics of MIMO-OWC are 
experimentally uncovered and our established criterion are verified. 

VI. Conclusion and Discussions 

In this paper, we have established a general criterion on the full-diversity space coded 
transmission of MIMO-OWC for the ML receiver, which is, to our best knowledge, the first 
design criterion for the full-diversity transmission of optical wireless communications with 
IM/DD over log-normal fading channels. Particularly for a 2 x 2 case, we have attained a 
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closed-form solution to the optimal linear FDSC design problem, proving that RC is the optimal 
among all the linear space codes. Our results indicate that the transmission design is indeed 
necessary and essential for significantly improving overall error performance for MIMO-OWC. 
However, the design criterion and the specific code constructions for MIMO-OWC presented in 
this paper are just initiative. Some significant issues are under consideration: 

1) Our proposed criterion can be applicable to any non-linear designs. It remains open 
whether there exists any better non-linear space code than RC. 

2) It has been shown in this paper that the space dimension alone is sufficient for full large- 
scale diversity. Here, a natural question is: what kind of benefit can be obtained if space- 
time block code designs are used for MIMO-OWCI 

3) Like MIMO techniques for RF communications, what is the diversity-multiplexing tradeoff 
for MIMO-OWC? 


Appendix 


A. Proof of Lemma 


Our main idea here is to split the whole integral in @ into two parts by properly choosing 


r > 0 such that its dominant term can be extracted. In other words, the average SEP [431 can 
be rewritten by 


PAp) = jfi(^\!H(h)dh 

+ Q(fF!k^fH{h)dh. (15) 

Now, we select r to satisfy that when p —)■ cx), r —)■ 0, to fragment (0, oo) into (0, r) and (r, cx)), 
adaptively with SNR. This fragmentation is to find the dominant term of (p) in the high SNR 
regimes by giving the upper-bound and the lower bound of Q fn {h) dh, and the upper- 

bounds of Q j fn (h) dh and then, examining their asymptotical behaviors related with 
r. Temporarily, we take this fragmentation for granted and then, will explain the essential reason 
later on. 

1) Upper-bound of SEP over (0,r).‘ We integrate the first part of SEP in ( [T5] ), which is 
denoted by Pr (p). Notice that when p —)■ cx, r —>■ 0, and in this case, fn (h) is monotonically 
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increasing over (0,r). Henee, Pr (p) is upper-bounded by 

/ 


Pr (p) < fHir) Q(^]dh 
Jo \ ^op 


< 


fnir) 


exp 


ph^ 


2P2 

op 


dh 


( 16 ) 


where the last inequality is obtained by using the Chernoff bound on the Gaussian tail integral. 
Furthermore, we ean upper-bound the last integral of ([T^ by 


exp 


'2P2 

op 


dh = PopV^p M o - Q 


P 


op 


< 


Pyp a/ _ 1 


-P 


(17) 


whieh produees 


p.(p) < 


Pop exp ^ 


p 2 exp 


(In r + cr^)' 
2cr2 


(18) 


2) Lower-bound of SEP: At the same time, when p —)■ cx), we have ^ <t. This inequality 
allows us to lower-bound P^ (p) by 

'y/ph' 


Prip) = 


Q 


P 


> 


> 


In p 


Q 


In p 


Q 


op 

y/ph 

P 

op 

y/pr 


fn {h) dh 


fn {h) dh 


fnih) dh 


^Poplnp 

where the last inequality follows from the monotonieally deereasing property of Q-funetion 
Then, by integrating fn {h) over fo, j, we have 


(19) 


Pr (p) > Q 


y/pr 

Pop In P 




— Inr + lnlnp\ 


Q 


> 


dp-r 

Pop In p / (T 


( 


2 ^ lnl^i£ 


where the last inequality is obtained by using 

1 


Q (x) > 


exp 


X 


a 


In^ 

In p 




2a2 


( 20 ) 


2y/^X 


V2 


exp — , X > —-. 


( 21 ) 
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3) Upper-bound of SEP over (r, cx)); We now turn to the seeond part of Pg (p) in ( [T5] ), whieh 
is denoted by Pr(p). Letting fn (h) = 0 produees the point ho satisfying \/h > 0, fn (h) < 

exp(^) 

fn (ho) = With this, Pr (p) can be upper-bounded by 

Pr{p) < fniho) Q(^\dh 


Lfsing Chernoff bound on the Gaussian tail integral gives us 


Pr{p) < ^fniho) / exp 


ph'^ 


2P2 

op 


dh 


= -fH{ho)V^P 


exp , „ 


p ' exp 


Pc 

pr" 


2P2 

op 


( 22 ) 


4) Determination ofr: Thus far, we have attained three bounds, respectively shown in ( [181 ), 
(|20j) and (|22|). Now, we examine their tightness related to r. We can see that ([T8]) and (|20|) have 


the same exponential term 




. To make the upper bounds in (181 to approach the upper-bound 


in ( |22| ) as tightly as possible, we seek to select such r that 

-2 (lnr + a2)' 


flL 

2P2 

op 


2a2 


(23) 


An explicit solution to ( [23] ) with respect to r is difficult to obtain. Accordingly, we propose to 
approximate the solution to ( [23] ) by 

Fop In p 


T = 


2Vp^ 


(24) 


which satisfies the following asymptotical equality. 


pP 


lim — "" 3 = lim 

p—^oo (Inr+tT^) p—^oo 

2a^ 


{InpY 




= 1 


Finally, combining ( [T8| ), ( |2()j ), ( |22| ), and ( |24| ) leads to the fact that there exist three positive 
constants Gil, Cnu and Ci 2 u independent of p such that 


G 


IL 


lnp-ln|r 

-‘op 


4.2 exp 


( 6np-lng^'^ 

8ct2 


< R(p) 
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<-(lnp) exp 


+Ci 2 uP 2 exp 


/ 

V 


li-i _i,-| 


4cr2 


8(7^ 


/ 


8(7^ 


where Cil = C'liu = \ and Ci 2 u = This is the eomplete proof of Lemma 


exp ^ 


(25) 


□ 


B. Proof of Theorem 

The eondition that any e G is unipolar without zero entry implies that if Vr > 0 and 
Ve G £, (e'^hj)^ < r^, then, we ean have hij < |^,i = 1,... ,M,j = 1,... ,N. After these 


preparations, we adopt the same teehniques as Subseetion III-A Temporally, assume that when 
p —)■ cxD, r —)■ 0. Similarly, P (s —)■ s) ean be adaptively fragmented with SNR as 

P(s^s)= [ P(s ^ s|H)/H(H)dH+ [ P(s ^ s|H)/H(H)dH. (26) 

The target in the ensuing subseetions is to give the asymptotieal bounds on P (s —s) following 
the similar procedures to the case of SISO-OWC. 

1) Upper-bound of PEP over (0, r); To begin with, let us process the first part of P (s —s) 
in ( |26l ) denoted by Pr (s —)■ s). We know when p —)■ cx), r —)■ 0. In this instance, fn (hij) 
is monotonically increasing over (^0, j and then, fn (hij) < Together with the 

Chemoff bound of Q-function, P^ (s —)■ s) can be upper-bounded by 


N M 


Pr (s ^ s) < - JJ JJ fn., 


T 


exp , I dH 


8NP^ 


.. (27a) 

j=li=l M-JI/ 7(erhp2<^2 

In addition, by Assumptionei, ■ ■ ■ , have the same signs. This result allows us to further 
upper-bound P^ (s —)■ s) by 


Pt (s s) < - 


n n 


i=l i=l 


T 


X 


exp 


PEf=iE£ie|/^., 


hij< 




dH 


(27b) 


Further, integrating the last term in (27b) produces 


N M 


Pr (s -)■ s) < - fn, 

j=l1=1 

N M 


N M 


nn^ 


Enn 

j=i *=i 


dP-f (_L) ( 
2 ■"‘«\\eU 


j = l i=l 
2 \ -| 




4Arp2 

op 


2 ^ V 4 Arp 2 


pP 


op 
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--MN 


o TT^ TT-^ ! 


P 


NP"^ 

^ op 


MN 

2 


exp 


N M /, 1 I I \ 2 

(lnr-ln|ej|) 


i=l i=l 


2a. 


(27c) 




2) Lower-bound of PEP over (0,r); To attain the lower bound of Pr (s —s), we need the 
following preparations. We observe that ee^ is rank-one with the only non-zero eigenvalue being 
Amax = e]. Then, hfee^hi < A^ax = I, ■ ■ ■, M. In other words, for all r > 0, 

it holds that 


N 


{H ; h, e = 1,... ,M} 

i=i 

C {H ; hj G , hfee^hj < r^, i = 1,..., M} 


Further, it is true that 


{H ; 0 < hij < 




-.,j = = 


N 


C{H;A^ax5^/iJ<r2,^ = l,...,M} 

i=i 

C {H ; hj G hf ee'^hj < r^, i = 1,..., M} 


(28a) 


(28b) 


Now, by (|28b[), we ean lower-bound Pr (s —s) by 


Pr (s —>■ s) ^ 


P(s ^ s|H)/h(H) dH 




> Q 


> Q 


\/Mp 


T 


\fWpT 


N ^2 / Jhij< 


/h (H) dH 


2NP^\Bpjj:li4.. 


N M 

nn« 

j=i 1=1 


In 




a. 


*2 


(29a) 


where the last inequality holds for high SNR sueh that In p > 1. Again, by pi] ), we arrive at 
the following lower-bound by 

\/MpT 


. exp (m 
Pr (S —)■ S) > —- ^ Q 


(dvr 


.MAT 




X • 


n N TT-/W’ 

,=i ni=i 




I^MN Inp-y/iVE^^iefc 


exp 


N M 

j=i i=i 


In 


r 


Inpt/iVELie: 


N „2 
k, 


(29b) 
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3) Upper-bound of PEP over (r, cxd); Now, we are in a position to analyze P (s —s) — 
Pr (s s), i.e., the seeond term of ( |2^ , which is denoted by Pr (s —)■ s). 

Similarly, for hij, ff (hij) = 0 gives the extreme point hijfl of fn (hij), i-e-, fn (hij) < 
fn {hijfl ), i = 1,..., M, j = 1,..., . Then, we have /h (H) < /h (Hq) and thus, Pr (s s) 
can be upper-bounded by 


Pr fs —s') < 


/h(Ho 


8NP^, 


l{eThif>T^ 


exp 


dH 


(30a) 


In fact, ( |28b[ ) also implies {H : hfee^h* > i = 1,..., M} C {H : htj > 




i = 


,M,j = 1,..., 'With this, Pr (s —>■ s) can be further upper-bounded by 


Pr (S —)■ S) < 


/h(Ho 


exp 






8NP^ 

Op 


dH 




N 2 
k=l 


i=l j=l 




4Arp2 

^ op 


Q 


pr^ 


4Ar2p2 
op 


exp 


< 


y^M y^JV 2 
2 ^j = \ ^ ij 


E N 2 

k=i 4 

NP^ 

op 




exp 


M pr"^ 

mp^ 

^ Op 


(30b) 


So far, we have attained three bounds and now, examine their tightness to select r in the 
following subsection. 


4) Selection of t: It is noticed that the exponential terms of ( |27c| ) and ( |29b[ ) are the same, 

2 sr^M -2 


i.e., (Inr)^ Ei=i J2j=i To mathch the bounds ( |27c| ) and ( |30b[ ), r is selected such that 


8A^P2 

^ op 


M N -2 
O', 


(top EE”" 


(31) 


i=l j = l 


Since a closed-form solution to ( [3T] ) is hard to attain, the solution to ( [3T] ) is approximated below 

(32) 


T = 


\ 


NP'^ 

^ op 

M 


M N 


2lnp 


i=l j=l 


For the selection in ( |32| ), the following asymptotical equality holds. 

MpP 

lim 




p^oo 


(lnr)^E.=iEL Y 


= lim 

p—^OO 


{In pY 


In -4 _In r y^ ^ ^ 

hPp I M 3^i=l 4^3=1 
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C^L (lnp)-“^^e 


_M N 

-MN __ ^*=1 ^J=i 


(lnp+ln(p|pn)-ln(MEf^ie|))" 


< P fs —s') 


Pi(s-S>s) 


<C2u (lnp)-^“^^e 


E lvl v-v i V 


r^ hlnf o') — In e^') o 

\ In^ p \ P / •?y ^-~yN In p 

+C^uP ^ e (33) 


P[ 7 l(s-s>s) 


■ V' 

Pu2{s^s) 


= 1 


Then, putting (27c), (29b), (30b) and ([32|) together, we are allowed to arrive at the fact that there 


exists three positive constants C 2 L, C 2 U and C^u, independent of p shown by ( [3^ , postponed 
to the top of the next page, where 


M N 




a. 


-2 


i=i j=\ 


(34a) 


n M y~\N 


(47r)"^^exp(-3fV) 


N 


C2L = UIE"' 


yk=l 


(34b) 


C 2 U — 


(NPl) 


MN 


r> TrAf TtA^ / 2 

2n.=i n,=i 


: exp 


a, .JNPip 

'Hm- 


(34c) 


exp 


C^u — 


2 


p2 

4^fc=iO 

NP'i 

op 


MN 

2 


(34d) 


2 n:Lin“i^u 

Now, we can see that in ( |3^ , P^ (s —)■ s) and Pu 2 (s —)■ s) have the same exponential term, 
exp ^ In^ p ), whereas the exponential term of Pui (s —)■ s) is exp ^ In^ ii)^) ’ decays 
slower than exp (—^ In^ p) against high SNRs. That being said, we have attained the dominant 
term. Pm (s —)■ s), of the upper-bound of P (s —)■ s). This completes the proof of Theorem □ 


C. Proof of Property^ 

For mi -L Pi, m 2 -L P 2 and m 3 ± P 3 , we have that 

+ P 3 . ^^2 ^ . 

- > — <;=> m2P3 — m^n2 > miP2 ~ m2Pi 

mi + m 3 m 2 
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This is indeed true, since we obtain min 2 — m 2 ni = 1 by Lemma and m 2 n 3 — m^n 2 > 1 from 
the assumption that ^ ns. 

In the same token, given mo -L no, mi ± ni and m 2 -L n 2 , we have that 


no + n2 ni 

- < — <;=> moni — mino > min2 — m2ni 

mo + m 2 mi 

This indeed holds, since we can attain min 2 — m 2 ni = 1 by Lemma and moni — mino > 1 
from the assumption that m < lU- xhis completes the proof of Property [T] □ 


D. Proof of Property 

First, Statement 1: — < < 111 . ^an be verified by 

“ m 1 m-i-Umo TTJ.o «' 


mi ' mi+m2 ' m 2 
ni + 02 02 


oil 


mi + m2 012 (mi + m2) \mi 

ni + 02 Oi 012 


(ni 02 


012 


< 0 , 


mi + m 2 mi (mi + m 2 ) \oi2 




mi 


When e ^e can have (mi/n - ni/ 12 ) > 0 and (m 2 /ii - 02 ^ 2 ) < 0. Hence, 


and 


and (m 2 /ii - 02/1 

(hi 

Oi+02\ 

V /12 

0 X 1 + 0 X 2 J 

(hi 

O 1 +O 2 \ 

V /22 

0 X 1 + 0 X 2 J 


(35a) 


(35b) 


Combining (35a) and (35b) leads us to the following facts: 


1) If e then, we have 

^ /i2 ’ /22 ymi’mi+m2y’ ’ 

-m 2 fii + 02 fi 2 > roifii - ni /12 > 0, - 0 x 2/21 + 02/22 > 0 x 1/21 - ni /22 > 0 


Thus F <c F 

Aiiuo, ± mini ^ rr 


holds. 


/ii hi r- i ni+n2 n2 


2) If 


/12 ’ /22 Vmi+m2’m2 


, — , then, we arrive at 

oxifii - ni /12 > - 0 x 2/11 + 02/12 > 0 
nTi /21 - ni /22 > - 0 x 2/21 + n 2/22 > 0. 


Hence, we have F^-^ni > F, 


m2n2‘ 


3^) If Ill = hi = ni-\-n2 ^ p 


/12 /22 mi+m2’ 

This completes the proof of Property 


m2n2- 


□ 
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E. proof of Property 

Vm satisfying m < lu < m ^hen e we have 

rr?,n * m.n J o rn.n mi m.o m.n ’ /n ‘ foo \ m.^ ' m,o / ’ 


Fmono, Fmsn^ > 0. Property [T] tells us that ^ and As 




m,o+m2 


a result, we ean have 


ni n 2 \ ^ / nQ+n 2 n 2 \ (^ ( no_ 2 l 2 _ 

mil m 2 J — ym,o+m'2’ m 2 J — ymo’ m 2 

ni n 2 \ ni+Uj \ ^ 

mi ’ m2 J — y mi ’ mi+ms J — y mi ’ m3 

Using Property]^ we ean attain F^ono > -^mjna and Fmmi < Fm^ns- The sueeessivity of ^ 
and ff ensures Fmim and F^ana are the last two minimum over (^, ^ )• This eompletes the 


m2n2 


proof of Property 


□ 


F. proof of Theorem 

Assume ^ and ^ are varying two sueeessive items of 5^2 p-i and denote x = jP and y = 
These notations produce /12 (1 + x) + /22 (1 + ?/) = 1. We firstly prove that if x,y E (^, , 

then the local solution to ( [TT] ) is determined by 

rii + 72,2 


_ /21 


mi ^ m 2 


X = y = 


fl2 — fxi — 


(36) 


mi + m 2 

mi + m 2 

2 (mi + m 2 + ni + 02 ) 

Letting (xmi - ni) {ymi - nf < (xma - 02) {ym2 - 02) yields x < InijlmfS^^"a-n'a) • 

In the same token, (xmi — ni) {ymi — ni) > (xm2 — 122) {y'm^2 — '^2) implies y < 

n2ixm2-n2)-ni{xmi-ni) and Property 1 ^ 

then ([TT]) can be equivalent to the following two sub-problems: 


. , t —Then, from Property |2|and Property |3[ we have that if x, y e 

m2[xm2—n2)—mi[xmi—ni) ’ ’ ^2 ’ 


1) If X < 


ni{ymi-ni)-n2{ym2-n2) 

mi{ymi—ni)—m2(ym2—n2) 


then ( [TT] ) is equivalent to the following problem 


max 

/lli/l2,/21,/22 


mini 


X < 


ni{ymi—ni)—n2{ym2—n2) 


S.t. < 


— mi{ymi—ni)—m2{ym2—n2)^ 


x,y E I ^ ^ 

’ ^ V mi ' m2 


(37) 


fll + fl2 + / 2 I + /22 = 1- 
Over the feasible region of x and with the increasing property of Fmim, we can have 

rii {ymi - ni) - n2 {ym2 - 02 ) 


F 


mini 


< mi 


mi {ymi — ni) — m 2 {ym2 — 72 - 2 ) 


- ^21 /12/22 {ymi - 72 i) 
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m2y {m2ni — min2) + n2 {min2 — m2ni) 
mi {ymi - rii) - m 2 {ym2 - n2) 

-f 12/22 {ym2 - 77 , 2 ) {ymi - ni) 


f 12/22 {ymi - ni) 


mi {ymi — rii) — m 2 {ym 2 — 77 , 2 ) 

where the last equality is attained by using min 2 — 7772771 = 1 and the equality holds 

This upper bound of is denoted by / (y) and 


when X = 

mi{ymi—ni)—m2iym2—n2) 

its derivative is 

/' (y) = 


((mi + 7772 ) 7 / - (771 + 772 )) 

(mi {ymi - ni) - m 2 {ym2 - 772 ))' 


((mi - m 2 ) 7/- (771 - 772 )) 


Letting /' ( 7 /) = 0 yields that y = 


ni—n2 


mi—m 2 


or 


ni+n2 

mi+m2 


. Now, we prove y = is the 

’ ^ ^ mi+m2 


is outside 

V mi ’ m 2 


mi—m 2 
77i 


We ean have 


-1 


maximum point. Firstly, we verify that 

771 - 77 2 _ 

mi — m 2 mi mi (mi — m 2 ) 

77i - 77 2 '^2 _ -1 

mi — m 2 m 2 m 2 (mi — m 2 ) 

From p^ , we know mi 7 ^ m 2 . Therefore, for the sign of mi — m 2 , we have the following 
two possibilities 

.1 ri-I rii-, n-I n ttt 1 .1 . / n 1 \ m / \ 

IS 


a) If mi > m 2 , then < ^ < We ean have that over ( f' (y) F 

^ mi—m 2 mi m 2 \ ^ 1 ’ mi+m2 J ^ 


positive and over 


/ ni+n2 712 


y miH-m2 ’ m 2 
ni—n2 


, ^ L /' ( 1 /) is negative. 


b) If mi < m 2 , then > — > —. Also, we can have that over . , , ,, 

mi—m 2 m 2 mi ’ y mi ’ mi+m2 / 


ni 711+712 \ 


/' ( 7 /) is positive and over 


, — ), /' (y) is negative. 

miH-m2 ’ m 2 / '' ® 


These observations indicate that y = is the maximum point. Combining y = 

^ TD.-X -t-m.o A O 


mi+m2 

and X = ”-lfa"^i~”-i)~”- 2 +m 2 -n 2 ) glygg ^ ^ y ^ ni+n 2 _ 
mi(^mi—nij—m2(ym2—^ 2 ) ® ^ mi+m2 




2) If 7/ < 


7l2{xm2—7l2)—7ll{xmi—7ll) 

m 2 {xm2—7i2)—mi {xmi —Tii) 


then ( [TT] ) is transformed into the following problem 


max 

flljl 2 j 2 lj 22 


m27l2 


y < 


7l2 {xm2—7l2)—7ll [xmi—Tli) 
m2{xm2—7i2)—mi{xmi—7ii) ’ 


S.t. < 


X, 7 / e I —, — 

’ ^ \ mi ’ m 2 


(38) 


/ll + /12 + / 2 I + /22 — 1- 


ni+n2 


In the same token, we can have the solution to (38) is x = 7/ = 

’ 1-r ^ mi+m 2 

Putting things together tells us that with the local solution to ( [TT] ) is 

/ll _ /^ _ 771 + 772 
/l2 /22 7771 + m 2 


(39) 
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It follows that Fmmi = Fm2n2 from Property 


Now, it is time to give the analytical FDSC. Combining with fu + /12 + /21 + /22 = 1 
produces /12 + /22 = ,Tii+^^2+ni+n2 • addition, using the geometrical and arithmetical inequality: 
^ for a, 6 > 0, we obtain 


F 


mini 


< 


f 12/22 {mix - ni) {miy - rii 
(/12 + f22f ( min2 - 


mi + m2 
2 




_ {min 2 - m 2 ni) 

4 (mi + m2 + ni + 72-2) 

where the equality holds if and only if /n = /21 = 


ni +112 


2(mi+m2+ni+n2) 

- Then, the local maximum value of the objective function is given by 


and fi2 = /22 = 


2(mi+m2+ni+n2) ‘ 


max min 6162 = -2- 

® 4 (mi + m2 + ni + 722) 

where we use min2 — m2ni = 1 from Proposition On the other hand. Proposition reveals 
mi + m2 + ni + 712 > 1 + 2 ^. This inequality implies 


maxmineie2 < - n- (41) 

F e ' " - 4(1 + 2 p)" 

where the equality holds when ^ ^ ^ or ^ ^ = 2 ^. Further, using fij = 1 

gives /ii = /21 = 1 ^. Finally, the global solution to ( [TT] ) is given by 


F™"^= 2 tW( 2 P l)' 

Until now, our focus is on the sequence before 1 / 1 . Using the reciprocity of the breakpoint 
sequence gives the final solution shown by ( [T?] ). This completes the proof of Theorem □ 
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